Effects of correlated disorder on the magnetism of double exchange systems 
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We study the effects of short-range correlated disorder arising from chemical dopants or local 
lattice distortions, on the ferromagnetism of 3d double exchange systems. For this, we integrate 
out the carriers and treat the resulting disordered spin Hamiltonian within local random phase 
approximation, whose reliability is shown by direct comparison with Monte Carlo simulations. We 
find large scale inhomogeneities in the charge, couplings and spin densities. Compared with the 
homogeneous case, we obtain larger Curie temperatures (Tc) and very small spin stiffnesses (D). 
As a result, large variations of ^2. measured in manganites may be explained by correlated disorder. 
This work also provides a microscopic model for Griffiths phases in double exchange systems. 



The interest for disordered magnetic systems 
such as thin magnetic films of transition metal 
alloys (Fe-Ni,Co-Ni,...), diluted magnetic semicon- 
ductors (Gai-^Mn^ASjGei-^Mn^,...), d° materials 
(Hf02,CaO,..) or manganites (Re x Ai_ x Mn03, where 
Re is a rare-earth ion and A an alkaline ion) has con- 
siderably increased during this last decade. One of the 
reasons is the potential of some of the materials to be 
incorporated in technological devices. Some of them play 
a very special role: systems which contain large scale 
inhomogeneities. Inhomogeneities can appear during 
the growth of the sample by molecular beam epitaxy 
for example but can also result from the interplay 
between many degrees of freedom (charge, spin, orbital, 
phonons). This is for example the case in manganites. It 
is known that manganites are strongly inhomogeneous 
at the nanometer scale: (i) large-scale structures in the 
charge density were seen by electron diffraction of thin 
films or tunneling spectroscopy 0]; (ii) evidence for 
inhomogeneous spin density was found in neutron diffuse 
scattering 0, or NMR 0; (iii) localized spin waves also 
suggest the presence of confining potentials [f|. There 
are also clear evidences of inhomogeneous structures 
above Tq, which were interpreted [a 0| as a Griffiths 
phase Q . Their microscopic origin is one of the central 
issues of the physics of manganites; it includes phase 
separation frustrated by long-range Coulomb interaction 



chemical disorder 



w long-range Coulomb lr 
HEK, polarons [HE1 



In this paper, we argue that the way the disorder is 
modelled is important to understand large-scale inho- 
mogeneous structures in 3d and to explain the Griffiths 
phase 0]. For this we consider a model where the disor- 
der is correlated at short distances. This model gives a 
possible explanation for the broad and multi-modal dis- 
tribution of NMR lines [411 , or the wide distribution of 
Curie temperatures Tq [3| and spin stiffnesses [IH mea- 
sured in different materials for the same carrier density. 
The microscopic origin of the correlated disorder could be 



chemical or polaronic. For instance, in Rex-aAsMnOa 
the dopant A 2+ which substitutes Re 3+ creates a strong 
Coulomb potential on its neighborhood and in particular 
on the eight nearest neighbor Mn sites surrounding it [lj| . 
This is the model of "color-centers" initially discussed by 
de Gennes (lfj| . Alternatively, local Jahn- Teller distor- 
tions can also be seen as a source of correlated disorder 
through "cooperative phonons", which can be mapped 
onto the same model. 

The 3d correlated disordered double exchange Hamil- 
tonian we consider reads, 



jh Si ■ & + i 



(i) 



where Uj = — t for nearest neighbors only, Si is a classical 
spin localized at site i (\Si \ = 1) and s*j = c\ a {S) a pCip, Jh 
is the Hund coupling which is set to be oo. The on-site 
potentials ej may correspond, in particular, to the chemi- 
cal substitution of Re 3+ by A 2+ defined by e, = en J2i x }> 
where the sum runs over the I nearest neighbour cations 
of the Mn site i (I = 1, . . . , 8), and en is the strength 
of the disorder. We choose randomly x cationic sites for 
A for which xj = 1 (otherwise x\ = 0). We emphasize 
that the disorder is correlated because one dopant af- 
fects simultaneously the 8 nearest neighbour Mn sites. 
With these definitions Ci takes the discrete values 0, eo, 
2eu,...,8eD- From stoechiometry, we would expect the 
hole density, to be equal to x, but in order to include 
the local Jahn- Teller distortion picture as well, we allow 
them to be different. 

The approach we use to study this model is in two 
steps. First, for a given configuration of disorder we di- 
agonalize Jl} in the real space, assuming a fully polar- 
ized ground-state at zero temperature. This allows to 
define an effective Heisenberg Hamiltonian for the clas- 



sical spins, H t 



ctr 



— zJ<ii> JijSi 



Sj, where the disor- 
dered couplings { Jij} are explicitly calculated in the limit 
Jh - » oo, using = (cKc Ji |)/2 [I?], 18|. In the 



2 



0.25 
0.2 



"0.15 



(a) 



— Mean-field, (effective model) 
Monte-Carlo (effective model) 

■ - RPA (effective model) 

■ Monte Carlo (MCI) 
O Monte Carlo (MC2) 
O Monte Carlo (MC3) 




J 




(A/tr 



FIG. 1: (Color online) (a) Tc as a function of hole density 
(clean case). Lines are obtained with the effective Heisen- 
berg Hamiltonian within Mean Field (dashed), RPA (dotted 
dashed) and Monte Carlo (continuous) treatments; and sym- 
bols from Monte Carlo simulations of the full double exchange 
model: MCI HH, MC2 [Ij and MC3 Q. (b) T c as a func- 
tion of the on-site potential width A for the Anderson dis- 
order (uncorrelated): The continuous line is obtained with 
Local RPA and symbols are from Monte Carlo [23| . 



second step, we diagonalize this Hamiltonian using the 
self-consistent local random phase approximation (SC- 
LRPA)[19]. It consists of decoupling higher-order spin- 
spin Green's functions in the equation of motion. This 
introduces the local magnetizations (Sf), which are self- 
consistent ly determined by using sum rules. Spatial fluc- 
tuations due to disorder are thus treated exactly by solv- 
ing the equations numerically in real space. This proce- 
dure was shown to be reliable to study dilute magnetic 
semiconductors where the couplings were calculated ab- 
initio 19]. SC-LRPA provides an analytical expression 
for Tq and allows us to study much larger systems than 
those used in Monte Carlo. 

In Fig. [H we test this method by comparing Tc with 
that of Monte Carlo simulations for both the clean sys- 
tem (a = 0) 0j2l|Jp, [13| and the system with An- 
derson disorder 23, 2m. In the later case, Ci are uncor- 



related variables uniformly distributed within [— -y,-^-]- 
In Fig. QJl (clean case), the lines are obtained by study- 
ing H c g within a simple mean- field theory, Tq = 2 J 
(dashed line), RPA, T c iPA = 1.32J [U (dotted dashed 
line) and Monte Carlo, T£ IC 



Tit^N~ where the kinetic energy (K) depends on 
nh- For the clean system, we remind that T PPA = 1.32 J 
is obtained analytically using T PPA = g(X^ q Wql) -1 ' 
where E(q) = zJ(l — 7(q)) is the magnon dispersion, z 
the coordination number, and 7(q) = ~ X)r. e 4qr ' [25j |. 
This expression actually gives a very good approxima- 
tion of Tc; the error compared to Monte-Carlo is 8%. 
Now, the comparison with Monte-Carlo simulations of 
the full double exchange model (symbols) shows that the 
difference is within 10%, so that the two step approach 
is quantitatively reliable. Similarly, when Anderson dis- 
order is added, we have found that SC-LRPA gives an 
excellent agreement with Monte Carlo data (Fig. [TJa), 
stressing that not only the thermal fluctuations are well 
treated but also the spatial fluctuations due to disorder. 

From now, we consider the model with correlated dis- 
order, as discussed above. In Fig. f2fi, f2b, [2fc, we have 
plotted the magnetic couplings Jjj in a given layer for 
various hole densities (rih — 0.1, 0.3 and 0.5, respec- 
tively). They are calculated for a fixed concentration of 
randomly distributed impurities (color centers), x = 0.3, 
and for a disorder strength, cd = 0.15 W (W = 12 t is 
the bandwidth), that is chosen to be compatible with ab 
initio calculations [l6j |. At low density (Fig.f2^i), the cou- 
plings are extremely inhomogeneous in space: we observe 
large clusters of strong couplings, embedded in regions of 
weak couplings. The distribution function of {2Jy } (not 
shown) is peaked at w — 0.003i but has a very long tail 
up to a cutoff of — 0M (the average is J = — 0.02i). The 
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FIG. 2: (Color online.) Top row: real space picture of the 
magnetic couplings 2Jy of the effective model (correlated dis- 
order), on one layer of the 12 3 cube. The dark (resp. white) 
regions correspond to large (resp. weak) couplings. Bottom 

(s z ) 

row: real space distribution of Xi — limT^T c -^ir on the same 
layer. From left to right, nn = 0.1 (a),(d), 0.3 (b),(e) and 0.5 
(c),(f). Parameters are x = 0.3 and to = 0.15VF. 

regions of strong couplings correspond to hole rich re- 
gions with metallic properties embedded in a hole poor 



3 



matrix which is expected to be insulating, thus leading 
to phase separation. This tendency will be reinforced 
if antiferromagnetic superexchange couplings are taken 
into account, the hole poor regions will become antifer- 
romagnetic or canted, as observed at very low dopings 
(droplets in a canted matrix) 0]. For Anderson disor- 
der, we do not have well-defined nanoscale regions in 
3d 24 1, unless cooperative phonons were included [27j ]. 
As the concentration of holes increases, the size of the 
regions of large couplings increases, and the system be- 
comes less inhomogeneous. In this respect, close to half 
filling (rih — 0.5), the nature of the disorder becomes 
less important, as we shall see. The reason is that car- 
riers with short Fermi wave length are less sensitive to 
the details of the disorder. Spatial inhomogeneities in 

the magnetization near Tc are directly seen in the dis- 

(S x ) 

tribution of A; = lim?^ j- c i— , where m is the averaged 
magnetization (Fig. \2§ . For a nearly homogeneous state, 
A, is close to 1, as seen in Fig. [2j. At low densities, we 
see a very inhomogeneous texture of Xi (Fig. [2ji), with 
local droplets with Ai as high as ps 4 — 5, surrounded by a 
region with very small local magnetizations. In this case, 
the distribution of the magnetizations is multi-modal. In 
between (Fig. [2b) , the droplet increases in size and A^ ps 2 
is reduced with respect to Fig. [2^, the distribution has 
only one broad peak. These results resemble NMR results 
where multimodal distributions occur at low dopings and 
get broader for higher doping 

In Fig. 02 we give Tc averaged over at least 100 dis- 
order configurations (symbols). To see clearly the role 
of the inhomogeneities, we have also indicated what Tc 
would be if we replace all couplings by their average, de- 
fined by J = J ij = (Jij)dis (lines). The results 
are almost identical for nh close to 0.5 but strongly dif- 
fer otherwise. Similarly we have found (not shown) that 
for Anderson disorder, Tc is also extremely close to that 
of the homogeneous system calculated with J. However, 
at lower hole densities where the couplings are strongly 
inhomogeneous (Fig. [2^), we observe that Tc is larger 
than that of the homogeneous sample. This happens 
because of the competition between large (percolating) 
clusters with couplings much stronger than J that tend 
to increase Tc and thermal fluctuations that reduce it. 
In particular, at rih = 0.1, Tc happens to be close to 
the mean- field result T C 4F = 2 J (dashed line), as the 
result of this competition. It is interesting to remark 
that this picture is different from the pure percolation 
picture where thermal fluctuations in the clusters wins 
and reduce Tc; the difference is that the distribution is 
much more inhomogeneous here. We note that our Tc 
are much smaller than that obtained in ref. [281 ] . where 
the same model was studied. The reason is that here 
both spatial and thermal fluctuations are treated beyond 
mean-field virtual crystal approximation. 

We now argue that this model gives grounds for a Grif- 
fiths phase [a, Q above Tc- As discussed in ref. 
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FIG. 3: (Color online.) Tc (symbols) as a function of hole 
density for the model with correlated disorder calculated by 
SC-LRPA (averaged over 100 configurations of disorder). Also 
given are the mean-field (dashed line) , RPA with all couplings 
identical Jij = J (continuous line). Parameters are x — 0.3 
and €d = 0.15W. Calculations are done for sizes 16 3 and 20 3 . 



correlations in the disorder should enhance the Griffiths 
phenomenon. Indeed, it is more likely to find large clus- 
ters with higher local "Curie" temperatures, as seen in 
Fig. [2^. We calculate this temperature Tc from the 
lowest eigenvalue e of Jy, using T G — + l)|e| 

0. Since this is a mean field estimation, Tc has to 
be compared to Tq F . For nh = 0.1, we find a large 
T G = O.llt ps 2.5T C IF . On the other hand, for n h ~ 0.5, 
the couplings are much more homogeneous, and we have 
found a much smaller region for the Griffiths phase with 
Tc = 0.13i ps T C AF . This is interesting because it shows 
that Tc is weakly sensitive to nh- Experimentally this 
phase seems to occur only in the structurally distorted 
phase at low dopings 0], which suggests that the origin 
of correlated disorder is the local Jahn- Teller distortions, 
a case that is also covered by the present model. In fact it 
is not clear from our study that we can exclude the chem- 
ical origin of the correlated disorder because the Griffiths 
phase shrinks as we increase the carrier density. A bet- 
ter treatment of thermal fluctuations could possibly lead 
to the complete disappearance of the Griffiths phase for 
larger dopings. 

We now discuss the effect of the inhomogeneities on 
the long wave length spin excitations at zero tempera- 
ture. Even in the presence of disorder, these excitations 
are well defined and characterized by a spin stiffness D 
that is calculated following [2^]. Experiments 
on various manganites show that the dimensionless ratio 
D/a 2 Tc (a is the lattice constant, taken to be 1 in the 
following) strongly varies with doping and takes values 
as small as 0.05 and up to 0.5 [15j. This is in contrast 
with the clean double exchange model, where D/T q is a 
constant equal to 0.755 (RPA for the s.c. lattice [25j), 
independent of the hole density. We argue that the mea- 
sured small values could be explained with the model of 
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occurence of a Griffiths phase observed for low dopings. 
In addition, we have found that correlations in the disor- 
der tend to increase Tq because of the presence of large 
clusters of strong couplings, and decrease the spin stiff- 
ness D. This results in very small ratios D/Tc, consistent 
with what has been measured in manganites but hardly 
compatible with Anderson disorder. The dimensionlcss 
ratio D/Tc appears as a good measure of the inhomoge- 
neous character of the magnetic state, a conclusion that 
may apply beyond manganese oxides. 

We thank M. Clusel, M. Hennion, P. Majumdar, Y. 
Motome, F. Moussa and S. Petit for stimulating discus- 
sions. O. C. thanks the ILL for hospitality. 



FIG. 4: (Color online.) Dimensionless ratio D/a 2 Tc of the 
spin stiffness to the Curie temperature as a function of the 
hole density rih, for the correlated and Anderson forms of dis- 
order. The width of the distribution of potentials was chosen 
to be the same in both cases: cd = 0.15W (a; = 0.3) and 
A = 0.80W. Experimental results flBI ] are divided by the 
lattice constant of LaMnOa squared (a — 3.9A). 



correlated disorder, but would require unrealistic large 
strength of the disorder in the uncorrelated case. Fig. [4] 
gives this ratio, calculated for both models of disorder. 
Note that to allow for a direct comparison, the width of 
the distributions of e$ is chosen to be the same. Close to 
rih = 0.5, D/Tc does not really depend on the model, 
reflecting the absence of the large scale inhomogeneities 
we discussed above. When rih decreases, however, the 
spin stiffness is dominated by large regions of weak cou- 
plings (Fig. 12k): at rih = 0.10, D/Tc is 3 times smaller 
than that obtained with Anderson disorder. In order to 
get such small values in the Anderson disordered case, 
one would need a value of A much larger than the band- 
width, which would be difficult to reconcile with ab ini- 
tio estimations 



16j | , o n one hand, and would tend to 



localize all carriers 24j on the other hand. In Fig. [|] we 
have also compared our calculations directly with avail- 
able experimental values [ll|. First we note that the 
overall quantitative agreement does not mean that ejj is 
quantitatively determined for manganites because other 
interactions have been neglected here (there is anyway 
a distribution of ratios for a same doping, see LaSr and 
LaPb in Fig. [U for instance which could be explained by 
different en)- Nevertheless, it is interesting to see that 
the trend as function of the hole density is already well 
captured by taking disorder into account, and that a rel- 
atively small amount of correlated disorder leads to very 
small values of D/Tc, contrary to what would be needed 
in the Anderson case. 

To conclude, we have found that short-range correlated 
disorder creates large scale spin and charge textures, par- 
ticularly inhomogeneous at low dopings. Our study sug- 
gests that describing the disorder in a more realistic man- 
ner may be a key point to understand experiments, as the 
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